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Abstract 



In this paper, we present a universal control technique, the non-holonomic 
control, which allows us to impose any arbitrarily prescribed unitary evo- 
lution to any quantum system through the alternate application of two 
well-chosen perturbations. 

1 Introduction 

Quantum control is very topical and interests many different fields of contem- 
porary physics and chemistry, such as Molecular Dynamics in laser fields and 
Quantum Optics 121 0] |S] . A few examples of control of the quantum state 
by conditional measurements 1 by adiabatic transport or by unitary evo- 
lution have been already proposed for the particular quantum system of atoms 
interacting with quantized electromagnetic field in a single-mode resonator. In 
parallel, a theoretical framework of quantum control has been built up : in par- 
ticular, several classes of problems have been distinguished such as the control of 
quantum evolution, quantum state, or density matrix ; besides, the mathemat- 
ical necessary conditions for their feasibility have been derived in the context 
of the theory of Lie groups (HHH. Moreover, theoretical methods, mainly based 
on optimization techniques, have been designed to exhibit the explicit values of 
the control parameters. 

Most of these traditional approaches rely either on a known or intuitively 
guessed particular solution which can be further optimized with respect to a 
given cost functional, through variational schemes ^Ol- By contrast, identify- 
ing the convergence domain of the standard algorithms through straightforward 
exploration of the entire multidimensional space of the control parameters im- 
plies numerical work, the complexity of which grows exponentially with the 
dimensionality of the system. Though, as we shall see in this paper, an original 
method, called non-holonomic bang-bang control UT], solves this problem for 
complex enough quantum systems so that no separable subsystems can be sin- 
gled out : in terms of classical mechanics, such systems do not have holonomic 
constraints. The physical idea is to alternately apply two distinct perturbations 
Pa and Pb during pulses, the timings of which play the role of control parameters 
and are determined by solving the " inverse Floquet problem" . Equivalently, one 
can fix the duration of the pulses and choose the strength of the perturbations 
as the free parameters one has to adjust in order to achieve the control objec- 
tive. Actually, the convergence of our algorithm results from an unsuspected 
simplification emerging from the Random Matrix Theory. Indeed, it relies on 
the algebraic properties of the TV*'' roots of the identity matrix, the spectra 
of which resemble to those of random unitary matrices which obey the Dyson 
distribution law. 

This paper is organized as follows. In the first part, we recall the broad 
context of the theory of quantum control. In the second part, we present the 
non-holonomic control technique, and provide all the algorithmic tools which 
allow one to implement it. 
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2 Quantum Control 



Quantum control is a very topical issue in contemporary physics. The needs for 
control are particularly sensitive in Molecular Dynamics and Quantum Optics, 
and are quite various : indeed, one should not speak of quantum control, but 
rather of quantum controls. Actually, four different types of problems have 
been identified in the literature ^lE]: the control of pure state, the control of 
density matrix, the control of observable and, finally, the control of the evolution 
operator. For each of these problems, the goal is the same, that is to impose 
the considered characteristics an arbitrarily chosen value. 

To achieve a control objective, one has to perturb the system, since its 
natural evolution usually results in too restrictive a dynamics. The control 
Hamiltonian H (t) comprises thus the unperturbed Hamiltonian Hq as well as 
M Hamiltonians of the form Cm (t) Pm , which can stand for the interaction 
Hamiltonians of the system with M classical fields, the amplitudes of which, 
denoted by Cm (t), are controlled by an external operator 



The functions {Cm (t)} play the role of the control parameters one has to adjust 
in order to achieve the desired control process. In other words, any problem of 
control can always be translated into the following form : for the physical system 
considered, perturbed through given fields of controllable amphtudes {Cm (i)}, 
one looks for the values of the control paramaters {Cm (t)} such that a specific 
characterics (quantum state, density matrix, observable, evolution operator) 
takes an arbitrarily prescribed value. 

All the objectives are not feasible. Some of them are forbidden by the uni- 
tarity of the evolution operator for closed systems : for example, one cannot 
change the eigenvalues of the density matrix through a Hamiltonian process of 
control. This kind of constraints is often referred to as kinematical constraints 
|13|. But there also exist dynamical constraints which stem from the algebraic 

properties of the Hamiltonians <^ Pm >■ Indeed, the evolution operator 



where T denotes the chronological product, belongs to the Lie group obtained 
by exponentiation of the Lie algebra generated by the operators <iHQ,iPm \- 



So, it appears clearly that the feasibility of a particular problem of control 
in a specific physical situation, defined by the Hamiltonians |-Pm|, is closely 
related to the properties of this algebra : for example, if one wants to completely 
control the evolution operator of a quantum system, that is if one wants to be 
able to give the operator U any prescribed value, one must perturb it in such a 
way that the operators |«i?o, iPm^ generate the whole Lie algebra u (N) which 
provides, through exponentiation, the whole Lie group U {N) jHl^j- Necessary 
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mathematical conditions also exist for the other types of control problems and 
can be found in literature |12| : these conditions are obviously weaker than the 
previous one, since the evolution controllability automatically implies all the 
other ones. 

The feasibility of a control problem can thus be decided through mathemat- 
ical criteria established in the context of the Lie group theory. But the explicit 
values of the control parameters achieving the desired control objective still 
remain to be found. In other words, once the existence of a solution has been 
proved, one has to find it explicitly. Different methods have been proposed, such 
as optimal control |15l I16L |2] , which mainly employ optimization techniques : 
the idea is to write down a functional the optimization of which gives the best 
compromise between all the constraints. A purely algebraic approach ^7], based 
on the decomposition of the arbitrary desired evolution on the Lie group, is also 
possible, but rapidly leads to intractable computations as the dimension of the 
state space increases. 

In the next section, we describe a universal control technique, called non- 
holonomic control jllj . which allows one to control the evolution operator of 
an arbitrary quantum system through the alternate applications of two well- 
chosen perturbations. Two equivalent sets of control parameters are possible : 
one can choose to freeze the amplitudes of the perturbations and let the pulse 
timings play the role of free parameters, or, on the contrary, freeze the pulse 
timings and adjust the amplitudes on each control pulse. We shall present the 
algorithm which allows us to compute the explicit values of control parameters, 
and emphasize the algebraic reasons which assure its convergence. 

3 Control of the evolution through non-holonomic 
control 

Let us consider an iV-dimensional quantum system of unperturbed Hamiltonian 
Hq. Our goal is to control its evolution operator U, i.e. to be able to achieve 
any arbitrary evolution U arbitrary S U [N). To this end, we alternately apply 
two physical perturbations, of Hamiltonians Pa and Pb, during N"^ pulses the 
timings of which are denoted by {r^ = — tk-i}k=i m (^o = and tjy^ = T 
correspond to the beginning and the end of the control sequence, respectively) 
: the total Hamiltonian takes the following pulsed shape (cf fig. ^1 

H{t) =Hn+Ca (t) Pa + Cb (t) Pb 

where 

Ca (t) = 1, Cb (t) = and H{t) = Ho + Pa ^ Ha for t E [t2k,t2k+i] 
Ca (t) = 0, Cb (t) = 1 and H{t) = Ho + Pb ^ Hb for t e [tsfe-i, i2fc] 
for fc = l,...,7V^ 
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Figure 1: Pulsed shape of the control Hamiltonian. 



and the total evolution operator is 



. e " 



where we have implicitly assumed that N is even. 

Our control problem can thus be translated into the following form : given 
U arbitrary S U (N) , an arbitrary unitary operator, we want to find a time vector 



V tn^ J 



such that 



U T ^ — U arbitrary ■ 



(1) 



As we said previously, for a solution to exist the operators |i-ffa, t-Hj^ must 

generate the whole Lie algebra u (N) : this is called the " bracket generation 
condition" . This property can be checked directly as long as the dimension N is 
not too big : one simply computes the commutators of all orders of iHa and iHb 
and stops as soon as they generate u{N). But when N becomes large, direct 
computation is intractable. In that case, one can simply check the following suf- 
ficient condition, suggested by V. Kac, according to which the system becomes 
non-holonomic, that is completely controllable, when the representative matrix 
of Hb in the eigenbasis of Ha has no off-diagonal zeros. 

Once the previous criterion is checked, one has to compute the time vector 
T solution of Eq. . The method consists first in determining the time vector 
T*(°) such that 

u = /, 

and then iteratively approaching the time vector through a Newton-like tech- 
nique. 



5 



The straightforward way to compute T*^'') would be to minimize the func- 
tional 

F(r) = lic/(y)-/d||2 

with respect to t . However, F presents many local minima which make its 
optimization uneasy. Nevertheless there exists an alternative method based on 
the algebraic properties of the TV*'' roots of the identity matrix. The idea is to 
look for TV parameters {Tk}k=i...N such that 



(2) 



where /« is an root of the identity matrix, i.e. a matrix the eigenvalues 



of which coincide with the N*^ roots of the unity 
words, /« is a matrix of the form 



in other 



/« = M- 



1 

e'T? 










•M, 



where M is a unitary matrix. To compute the T/t's, we use the following alge- 



braic property : if Pjj{X) = X^jLo denotes the characteristic polynomial of 
a unitary matrix U, then X^^Lo — ^ and the equality is achieved iff U is an 
TV*'* root of the identity matrix, up to a global phase factor. To obtain the T/c's, 



one thus computes the characteristic polynomial P(A) = ^ 
of the matrix product 



N 



{{Tk}k=i...N)Xi 



and minimizes the function = X]j=o I'^i ({'A;}fe=i...jv) | to 2 with respect 
to the Tfc's. This minimization turns to be quite easy, due to the fact that 
a generic unitary matrix is very close to an A^*'* root of the identity. In fact, 
numerical work shows that in about 30% cases of randomly chosen timings {Tk} 
the standard steepest descent algorithm immediatly finds the global minimum 
Fn — 2. This fact has roots in the Random Matrix Theory. Indeed, according 
to Dyson's law, the eigenvalues of random unitary matrices tend to repel each 
other, and are thus very likely to be almost regularly distributed on the unit 
circle, as those of an N*'^ root of the identity, as shown in Fig[21 In other 
words, in the space of N x N unitary matrices, the /« matrices are present in 
abundance, and can be reached from randomly chosen point by small variation 
of the timings. 

Finally, we define the time vector l^^") corresponding to the identity matrix 
by simple repetition of {Tk} 



\+U-l)N 



T, for 



,TV, 



(3) 
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Figure 2: Spectrum of a random unitary matrix (white circles) resulting from 
the repulsion of the eigenvalues on a unit circle is shown vs the eigenvectors 
of -/V*'' root of the identity matrix (black circles) multiplied by an unimportant 
phase factor. 



and checks that indeed 



In In 



N times 



up to an unimportant global phase factor. 

We now have to iteratively determine the time vector V from V^^^. Let us 
first consider the case of a target evolution close to the identity : in that case, 
Uarbitrary Can be Written under the form 

Uarbitrary = Ue = exp(-iHe), (4) 

where W is an 8 x 8 bounded (||7i|| < 1) dimensionless Hermitian Hamiltonian, 
and e > a small parameter. We then look for the variations St^, determined 
to first order in e by the linear equations 




Once has been calculated through standard techniques of linear algebra, we 
replace T*'"' by V-*^^ + (5V and repeat the same operation until we obtain V 
which checks J7 ( V) = Uarbitrary at the desired accuracy. 
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If the evolution UarUtrary — Ue IS uot close to the identity, that is if e 
is not small, one has to divide the work into elementary paths on which the 
previous method converges. To this end, we consider an integer n > 2 such that 

(^Uarbitrary^ = [/jL is attainable from / through our iterative algorithm, and 
determine in this way the associated time vector ^(i) which checks 

1 

Taking (^Uartitrary^ = U-^ as our new target, we repeat the same algorithm 
to compute t such that 

and so on. We progress in this way as long as our algorithm converges : in 
general, it stops at a value n* > 1, for which the system Eq|(5l) has no solu- 
tion. Then, we keep the time vector ~t ^^'^ and simply repeat the same control 

sequence n* times to achieve the desired evolution 



U 



- U, 



arbitrary • 



To conclude this section, we provide an equivalent form of our method. 
Indeed, in the previous paragraphs we fixed the amplitudes of the perturbations 
once for all and considered the pulse timings as our free control parameters. 
But we also might have chosen to apply iV^ pulses of same duration t = 
where T is the total control sequence duration, and taken the amplitudes as 
our free control variables. In other words, we might have applied the following 
Hamiltonian 

H2k+l = Hq + C2k+lPa 
H2k = Ho + C2kPb 



H{t) = 



(6) 



where the amplitudes {Ck}^^^ ^2 play the role of adjustable control param- 
eters. The evolution operator would have thus taken the form 



U{{C,}) 



. . e 



(7) 



and our problem of control would have boiled down to finding the vector C 
( ^1 \ 



such that U 



(^) 



arbitrary- 
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The method remains the same as previously. First, one determines the 
parameters {ckj^^i jv such that 

e--l(^o+^JvP!,)r . ^-^{Ho+CN-iPa)T _ _ _ g" I (ffo+ci Pa )r ^ ji^ 

by minimizing the functional X]^o I'^i {{ck}k=i...N)f to 2, where {aj} denote 
the coefficients of the characteristic polynomial of the matrix product 

and we set 

Cl%_^j^ = ci for i,j = l,...,N, (8) 

such that U = /. 

Then, one iterativcly approaches the desired vector C in the same way as 
previously. For a target evolution UarUtrary = 11^= exp(— iHe) close to the 
identity {i.e. for e small), one computes the variations SCk to first order in e by 
solving the linear equations 



Then we replace C by C + 5 C and repeat the same operation, and so on, 
until we get C which checks U (^C^ = UarUtrary with the desired accuracy. 
For a target evolution U arbitrary = far from the identity {i.e. for e fi- 

nite), one considers an arbitrary integer n > 1 such that (Uarbitrary^ " ~ 
can be reached from the identity through the previous iterative algorithm, and 
computes C^"^ which checks 

Then, one calculates in the same way the vector C such that 

p(^))=f/^, 



and so on, until one obtains the limiting value n* > 1 beyond which the algo- 
rithm fails to converge. Finally, one gets the desired evolution by repeating the 
same control sequence n* times, the amplitudes of which are given by the vector 



4 Conclusion 

In this paper, we presented a universal method for controlling the evolution 
operator of any quantum system. This control scheme, called the non-holonomic 
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control, consists in alternately applying two physical perturbations which check 
the bracket generation condition. The role of tunable control parameters may be 
played either by the durations of the interaction pulses or by the strengths of the 
perturbations during these pulses. The calculation of the control parameters is 
achieved through an algorithm, that we presented here, the convergence of which 
is assured by remarkable properties of the iV*'' roots of the identity matrix. 

In the following paper, we show how the non-holonomic control technique 
can be used in the context of quantum computation in order to build controlled 
quantum devices. 
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